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Abstract 

Trapping of two-armed (m = 2) vertical p-mode oscillations in relativistic disks are 
examined. The disks are assumed to be isothermal in the vertical direction, but are 
truncated at a certain height by the presence of corona. The same issues have been 
examined in a previous paper (Kato 2012a). In this paper, unlike the previous paper, 
however, we do not use the approximation that the oscillations are nearly vertical, 
but limit to a simpler case of no magnetic field. As in the previous paper, the results 
suggest that the two basic oscillation modes [both are the fundamental (n = 1) in 
the vertical direction but in the horizontal direction one is the fundamental (n r = 0) 
and the other the first overtone (n r = 1)] correspond to the twin kHz QPOs. Second, 
the oscillation mode which is the first overtone (n = 2) in the vertical direction and 
the fundamental in the horizontal direction (n T = 0) will correspond to the horizontal 
branch oscillation (HBO) of Z-sources. The latter suggests that the horizontal branch 
of Z-sources is a sequence of temperature change in disks whose vertical thickness is 
strongly terminated. The temperature increases leftward along the sequence from the 
apex between normal and horizotal branches. 

Key words: accretion, accrection disks — horizontal branch oscillation — kHz 
quasi-periodic oscillations — neutron stars — two-armed disk oscillations — X-rays; 
stars 

1. Introduction 

KHz quasi-periodic oscillations (kHz QPOs) observed in neutron-star low-mass X-ray 
binaries (NS LMXBs) will give an important clue to understand the innermost structure of 
relativistic disks of NS LMXBs, if their origin is clarified. So far, many models of kHz QPOs 
are proposed. Among them trapped disk oscillations are one of promising candidates, and 
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trapping of p-, g-, and c- mode oscillations in relativistic disks has been extensively studied 
(see Wagoner 1999, Kato 2001, and Kato et al. 2008 for review, and see Fu & Lai 2009 and 
Lai & Tiang 2009 for related topics). 

In addition to the above trapped oscillations, Kato (2010, 2011a,b, 2012a,b) suggested 
that trapped, two-armed (m = 2), vertical p-mode oscillations in geometrically thin (magne- 
tized) disks are one of possible origins of kHz QPOs observed in neutron-star X-ray sources. 
In this model two trapped oscillations, both of which are fundamental modes (n = 1) in the 
vertical direction but one is the fundamental (n r = 0) and the other the first overtone (n r = 1) 
in the radial direction, are regarded to be the twin kHz QPOs. Frequencies of these two os- 
cillations change with correlation by time changes of disk temperature, disk vertical thickness 
and toroidal magnetic fields in disks. This correlated change of the frequencies is considered to 
be the cause of the observed correlated change of the twin kHz QPOs (Kato 2011b, 2012b). 

In calculating the frequencies of trapped, two-armed p-mode oscillations, we have as- 
sumed that the disks are isothermal in the vertical direction and have toroidal magnetic fields. 
Furthermore, for simplicity, the toroidal magnetic fields are assumed to be distributed in the 
vertical direction so that the Alfven speed, ca, is constant in the vertical direction (Kato 2011a, 
2012a). Even in this simplified case, semi-analytical calculations of frequencies of trapped ver- 
tical p-mode oscillation are still complicated. Hence, in previous studies (Kato 2011a, 2012a) 
we have assumed that the oscillations are nearly vertical and the motions in the radial direc- 
tion associated with the oscillations can be treated as small perturbations over the vertical 
motions. The approximation is allowed in qualitative arguments, but not accurate enough to 
do quantitative comparison of our model with observations. 

Hence, in this paper, returning to a simpler case, i.e., to the case of non-magnetized 
disks, we examine the trapping of two-armed vertical p-mode oscillations, since in this simpler 
case the frequencies of trapped oscillations can be calculated analytically without introducing 
the approximation that the oscillations are nearly vertical. The results show that even in this 
case of non-magnetized disks, the observed correlation of twin kHz QPOs and the correlation 
between kHz QPOs and horizontal branch oscillation can be both described by moderate time 
variations of disk temperature, if the disks are strongly truncated in the vertical thickness. 

2. Disk Models and Parameters Describing the Disks 

Distinct from previous papers by Kato (2011a,b, 2012a,b), no magnetic field is assumed 
in this paper. Except for this, our disk models are the same as those adopted in previous papers. 
That is, the unperturbed axisymmetric disks are vertically isothermal. Thus, the vertical 
hydrostatic equilibrium gives that in the vertical direction the density po(r,z) is distributed as 



where the scale height H(r) is related to the isothermal accoustic speed c s (r) and the vertical 
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epicyclic frequency fl±(r) by 
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Here, r is the radial coordinate of the cylindrical coordinates (r,(p,z) whose 2-axis is perpen- 
dicular to the unperturbed disk plane and the origin is at the disk center. It is noted that the 
Newtonian formulation is adopted in this paper, except when the radial distributions of fl(r), 
n_i_(r), and n{r) are considered, where n{r) is the radial epicyclic frequency and Q(r) is the 
angular velocity of disk rotation and taken to be the Keplerian angular velocity, Q^(r), since 
geometrically thin disks are considered. 

We assume that the isothermal disks are terminated at a certain height and surrounded 
by hot, low-density corona. The height of termination, z s (r), is taken to be a free parameter. 
Hereafter, we introduce a dimensionless parameter, r) s , defined by t] s = z s /H, which is a function 
of r in general, but taken to be a constant independent of r, since the trapped region of 
oscillations is not wide, except for the case of small rj s , as is shown below (see figures 2 and 3). 

The disk temperature is also a parameter describing disks. As a reference, we adopt the 
temperature in the standard disks where the gas pressure dominates over the radiation pressure 
and the opacity comes from the free-free processes. In this case, if the conventional viscosity 
parameter a is 0.1 and the mass accretion rate normalized by the Eddington critical accretion 
rate is 0.3, the square of the isothermal accoustic speed, (c 2 (r)) , is given by (e.g., Kato et al. 
2008) 

/ M \ -V5 / r \ -9/io 

(c 2 (r)) = 1.79xl0 16 (-^ pi cm 2 s- 2 , (3) 
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where M is the mass of the central source and r g is the Schwarzschild radius defined by r g = 
2GM/c 2 . In this paper we adopt j3[= c 2 /(c 2 ) ] as a dimensionless parameter describing the 
temperature distribution in disks. (3 is, in general, a function of r, but taken to be constant 
independent of r, since the trapped region of oscillations is generally not wide, as mentioned 
before. 

In summary, we specify our disk models by two parameters: 

2 

and two additional parameters M/M & and a*, where a* is the spin parameter of the central 
source (0 < a* < 1.0). 

3. Trapped Oscillations and Their Frequencies 

Wave equations describing isothermal oscillations in vertically isothermal and geometri- 
cally thin disks have been well studied (e.g., Kato 2001 for review). It is convenient to introduce 
the variable hi defined by hi =pi/po, where p\ is the pressure variation over the unperturbed 
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pressure p and p is the unperturbed density. Then, under the approximation that the radial 
wavelength of the perturbations is smaller than the scale length of the radial variations of the 
equilibrium properties of the disk, we have a wave equation in terms of hi as (e.g., Kato 2001) 
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hi = 0, 



(5) 



(u — mVl) 2 — k 2 dr 

where the perturbations are taken to be proportional to exp[i(ut — mtp)], oj and m are being, 
respectively, the frequency and the azimuthal wavenumber of perturbations. In this paper we 
consider only the case of m — 2, but m is retained hereafter without specifying m = 2 so that 
the origin of m can be traced back. Let us introduce a dimensionless variable rj defined by 
rj = z/H. Then, equation (jSJ) is written as 
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Wave equations similar to or more general than equation §5§ have been solved by the 
WKB method, approximately separating hi(r,r)) as hi(r,rj) = g(rj)f{r,rj) (e.g., Silbergleit et al 
2001, Ortega-Rodrigues et al. 2008, Kato 2012a). In the present case, h\ is simply separated 
as hi(r,rj) = g{rj)f(r), and we have a set of two equations (e.g., Okazaki et al. 1987): 
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where K is the separation constant and can be obtained by solving equation ((Tj) with relevant 
boundary conditions, which is given in the next subsection. 

3.1. Oscillations in the Vertical Direction and Their Eig en-Frequencies 

In the case where the isothermal disks extend infinity in the vertical direction (i.e., 
r] s = oo), the boundary condition of finite wave-energy density at r] = ±oo requires that K is 
zero or a positive integer, say n, and 17(77) is a Hermite polynomial, ^(77) (Okazaki et al.1987). 
The integer n thus represents the number of nodes of hi in the vertical direction. In the present 
problem, we are interested in vertical p-mode oscillations and thus n starts from n — 1, not 
from n = (71 = corresponds to the horizontal p-mode oscillations and the oscillations of n = 
is outside of our interest in this paper.) Hereafter, we use n(= 1,2,3...) to denote node number 
of hi in the vertical direction. 

In the case of vertically terminated disks (?7 S 7^ 00), we impose boundary conditions that 
the Lagrangian pressure variation vanishes at 77 = ±r/ s , i.e., 5pi = at 77 = ±77 s , since in the 
coronae the perturbations propagate away rapidly because of high temperature. If we adopt 
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Table 1. Eigenvalue K„ tB of Vertical Oscillations 
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the approximation that the oscillations are nearly vertical there, this boundary condition can 
be reduced to (Kato 2012a) 

_£ = at 77 = ±r7 s , (9) 

where u z is the z-component of velocity associated with the oscillations. Since the ^-component 
of equation of motion is related to hi by i(u — mfl)u z = —dhi/dz, this boundary condition can 
be expressed in terms of g as 

dr] 

The eigenvalue K obtained by solving equation ([7j) with this boundary condition is expressed 
hereafter by K ns , since it depends upon 77 s and the number of node, n. The tilde is attached to 
K in order to avoid confusion with K n)S in Kato (2012a). The results are summarized in table 
1 for some values of n and 77 s . The eigenvalue K n s is related to K n>a obtained by the boundary 
condition dg/drj = by K n s = 1 + -K^sll] 



at ?7 = ±77 s . (10) 



3.2. Frequencies of Trapped Oscillations 

We solve equation ([8]) by the WKB method with K U}S given above. By introducing r(r) 
defined by, as Silbergleit et al. (2001) did, 

r(r)= [{[u-mn^f-K^r'^Xdr 1 , (11) 



where r ; is the radius where inner boundary condition is imposed, we can reduce equation 
in the form 

d 2 f 



dr 2 



+ Qf = 0, (12) 



where 



r]—+Kg = 0, (*) 



Let us denote the eigenvalue and cigenfinction of equation 

d 2 g dg_ 

drf dr] 

K niS and g u , respectively, when the equation is solved by the boundary condition dg/dr) = 0. Then, K ntS 
and g u are related to K n s and g given in this paper by K n ^ s = K n s — 1 and g u = dg/rj, respectively. This 
can be shown by substituting g u given above into equation (*) and using the fact that g is a solution of 
equation (*). 
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Fig. 1. Frequency - temperature relation of trapped two-armed (m = 2) vertical p-mode oscillations in 
isothermal disks. Two cases are shown where the disks extend infinitely in the vertical direction {z s /H = oo) 
and are terminated at the height of 2H, i.e., 77 s = z s /H = 2.0. Four oscillation modes are shown, i.e., the 
set of (n, n r ) arc (1,0), (1,1), (2,0) and (2,1). The thin curves are for the cases where the frequencies are 
calculated with the approximation that the horizontal motions associated with the oscillations are small 
perturbations over the vertical ones. In oscillations of n = 1, the frequencies are independent of 77 s . The 
mass M and the spin parameter a* adopted are M = IAMq and a* = 0. 

As r\ we take the radius of the marginally stable circular orbit and assume that h\ vanishes 
there, i.e, / = 0§| The other boundary condition is that outside of the capture radius, r c , which 
is defined by Q(r c ) = 0, the oscillations are spatially damped. Then, as the trapping condition, 
we have (see Silbergleit et al. 2001; Kato 2012a) 

) 1/2 ^Wr= / ^ 1 — 1 n ' s ± dr = vr(n r + 3/4), (14) 



/ Q 1 '\r)dr = 



C s 

where n r (= 0, 1,2, ...) is the node number of / in the radial region of n < r < r s . This equation 
gives the frequencies of trapped oscillations (see also the appendix), as functions of n, n T , (3, 
and ?7 S for given M and a*. 

4. Numerical Results on Trapped Oscillations 

Numerical results obtained by solving equation (fl~4|) are shown here mainly in the case of 
M = 1.4M with some values of a*, /3, and r] s . Among many oscillation modes of n(= 1,2,3,...) 
and n T (— 0,1,2,...), our main interests here are on oscillations modes of (n,n T ) = (1,0), (1,1), 
and (2,0), since we think that the former two will correspond to the set of twin kHz QPOs, and 
the third one will to the horizontal branch oscillation (HBO). 



The present boundary condition, / = 0, is consistent with the boundary condition adopted by Kato (2012a), 
which is u z = at r = r v This is because the z-componcnt of equation of motion is i(ui — mD,)u z = —dh\/dz 
and since dh\/dz = when h\ is separated as hi = g{rj)f{r) and / = 0. 
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Fig. 2. Propagation diagram for two oscillations of (n, n r ) = (1,0) and (1,1)- The horizontal lines show 
the radial range where oscillations are trapped. The thin lines are for the cases where the approximation 
of nerly vertical oscillations is adopted to calculate frequencies. f3 = Cg/(Cg)o = 10.0 has been adopted. As 
in figure 1, M — 1AM® and a* = are taken. 

Figure 1 shows the /3-dependence of frequencies of the trapped oscillations. Four oscil- 
lation modes of (n,n r ) = (1,0), (1,1), (2,0), and (2,1) are shown for two cases of ?7 S = oo and 
?7 S = 2.0 with a* = 0. For comparison, the results obtained by the approximation method that 
the horizontal motions associated with the oscillations are small perturbations over the vertical 
ones (Kato 2010, 2011a) are shown by thin curves. It is noted that in the case of n = 1, the 
frequencies are independent of rj s , since K n>s is 1.0 and free from r] s . Frequencies of oscillations 
decrease in the order of (n,n T ) = (1,0), (1,1), (2,0), and (2,1), as long as (3 < 10.0, and decrease 
as (3 increases. It is also noted that the frequencies derived by using the approximation of nearly 
vertical oscillations are systematically higher than those derived without the approximation, 
especially in the case of n = 1 . 

Figure 2 is the propagation diagram for two oscillations of (n,n r ) = (1,0) and (1,1) in 
the case of (3 = 10.0. Disk parameters adopted in figure 2 are the same as those in figure 1. 
For comparison, the oscillations derived by the approximation of nearly vertical oscillations are 
shown by thin curves. 

To understand the propagation diagram of figure 2, let us examine propagation and 
evanescent regions of an oscillation with frequency ujq. Since the propagation regions are spec- 
ified by Q > [see equation (Jl2|) ]. we know that the regions are (oj q — mVt) 2 — K n>s Vt\ > and 
(uo —mOt) 2 — k 2 < in the case where K n s > 1. The former is the propagation region of vertical p- 
mode oscillations and separated into two regions of loq — mfl < —KlJ 2 Vl± and ujq — mVL > K^ 2 ^l±, 
and the latter is the propagation region of g-mode oscillations and — k < uoq — mQ < k. In the 
case of vertical p-mode oscillations with n = 1, we have K n ^ s = 1. Hence, for oscillations of m = 2, 
ujq < 2Q — Q± is a propagation region, which is the region of our interest in this paper. The 
radius where uq = 2Q — Qj_ is the capture radius, r c . Outside of r c , there is an evanescent region 
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Fig. 3. Propagation diagram for oscillations of (n, n r )=(2,0) for three cases of rj s = oo, 2.0, and 1.1. Other 
parameters are the same as figure 2. 

which extends till the radius of ujq = 2Q — k (inner Lindblad radius). The region specified by 
2Q — k < u < 2Q + k is again a propagation region, but the vertical p-mode oscillations which 
are nearly trapped in the region of cuo < 2Q — Q± will penetrate only weakly to this propagation 
region. It is noted that inside the region of 2Q — n < u Q < 2Q + n there is the corotation point, 
ujo = 20, where the oscillations will be damped. In the present oscillations, however, their 
damping at the corotation resonance will be negligible, since wave leakage to the point will 
be negligible by spatial damping in the evanescent region, as mentioned above. The trapping 
condition Q14p determines the frequency and the capture radius which are shown in figure 2 
by thick horizontal lines. The thin horizontal lines show the case where the oscillations are 
approximated to be nearly vertical ones. 

Figure 3 shows the propagation diagram of oscillations of n = 2 and n r = in the case of 
(3 = 10.0. Three cases of r] s (i.e., r] s = oo, 2.0, and 1.3) are shown with other parameters being 
the same as those in figures 1 and 2. A purpose of showing this figure is to demonstrate that in 
the case of n = 2 the curve specifying the capture radius, i.e., the radius of uio = 2VL — if^/ 2 fij_ 
depends on 7/ s since K nfi depends on r] s . It is notices that when 7/ s is close to 1, the trapped 
region is rather wide, since K n ^ s is large and close to 4.0. 

Figure 4 shows the frequency- spin relation for four oscillation modes of (n,n v ) = (1,0), 
(1,1), (2,0), and (2,1) in the case of (5 = 10.0 with M = 1.4M . Two cases of rj B = oo and r/ s = 2.0 
are shown. As a* increases the frequencies of trapped oscillations increase, since the oscillations 
are more localized in the inner region as a* increases. 

5. Frequency Correlation and Comparison with Oscillations 

We take the standpoint that two oscillations of (n,n r ) = (1,0) and (1,1) with m = 2 are, 
respectively, the upper and lower kHz QPOs. In addition, we suppose that the (2,0) oscillation 
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Fig. 4. Frequency - spin relation for four oscillation modes of (n, n r ) = (1,0), (1,1), (2,0), and (2,1). In 
the case of n = 2, the relation depends on ?7 S . Two cases of rj s = oo and 2.0 are shown for the oscillation 
of n = 2. M and /3 adopted are M = 1.4M Q and /3 = 3.0. 

corresponds to the horizontal branch oscillation (HBO). These possibilities are examined here 
by comparing correlation curves derived by the model with observations. 

5.1. Frequency Correlation between (1,0) and (1,1) Oscillations 

We examine how the set of (1,0) and (1,1) oscillations move on a frequency-frequency 
diagram when f3 is changed from 1.0 to 10.0. First, two cases of a* = and 0.2 are shown in 
figure 5 with M = 1.4M Q . These correlation curves are superposed on the diagram showing the 
observed frequency-frequency plots of the twin kHz QPOs. It should be noted here that the 
correlation curves are free from ?7 S , since the frequencies of the n = 1 oscillations are independent 
of rj s . An increase of a* shifts the correlation curve in the upper-right direction. Let us next 
consider effects of mass of the central source on the correlation curve. Figure 6 is for M= l.8M@ 
and three cases of a* = 0, 0.2, and 0.4 are shown. This figure shows that an increase of 
mass shifts the correlation curve in the lower-left direction on the diagram, compared with the 
corresponding cases of M = 1.4M Q . 

From figures 5 and 6 we see that the observed correlation between the twin kHz QPOs 
of typical Z-sources can be well described by a* ~ if the mass of the central source is slightly 
larger than 1.4M Q . A larger a*, say a* ~ 0.2, is required if the mass is larger, say 1.8M . 

In previous papers (Kato 2011a,b, 2012a,b), we have calculated frequencies of trapped 
vertical p-mode oscillations and correlations among them in the case where the disks have 
toroidal magnetic fields. In these papers, however, we adopted the approximation that the 
oscillations are nearly vertical in the sense that the horizontal motions associated with the 
oscillations can be treated as small perturbations on the vertical motions. This approximation 
is qualitatively acceptable, but not always relevant to quantitative arguments. In order to 
see this more in detail, we show here the results derived by using the approximation in the 
particular case of no magnetic fields, since in this case we have derived the frequencies without 
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Fig. 5. The frequency correlation between (n, n T )= (1,0) and (1,1) oscillations for two cases of a* = and 
0.2. The value of (3 is changed from 1.0 to 10.0. M = 1.4M© is adopted. The plot of observational data 
of some typical Z-sources, taken from the figure of Abramowicz (2005), are overlapped on this figure. The 
straight line labelled by 3 : 2 is the line on which frequency ratio is 3 : 2. (Color Online) 




Fig. 6. The same as figure 5, except that M = 1.8Mq and three cases of a* = 0, 0.2, and 0.4 are considered. 
(Color Online) 

the approximation. The correlation curves of the (1,0) and (1,1) oscillations derived by using 
the approximation are shown in figure 7, when f3 is changed from 1.0 to 10.0 in two cases of 
a* = and 0.2 with M = \AM®. This figure should be compared with figure 5 obtained without 
the approximation. The comparison shows that in the approximate case the correlation curve 
shifts in the upper-right direction almost along the curve, compared with the case without 
the approximation. In other words, if an observed frequency correlation is described by the 
results obtained by using the approximation, we have larger values of /3 compared with those 
derived without the approximation. This can be understood from figure 1, where the frequencies 
calculated by the approximation are shown to be systematically higher than those calculated 
without using the approximation. Comparison of figure 7 with figure 5 also shows that the 
approximation gives systematically a larger a*, if other parameter values are fixed to be the 
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Fig. 7. The same as figure 5 except that the frequencies are calculated by using the approximation that 
the horizontal motions associated with the oscillations arc small perturbations over the vertical motions. 
(Color Online) 

same in both cases. 

5.2. Frequency Correlation between (1,0) and (2,0) Oscillations 

The frequency correlation between two oscillations of (n,n T ) = (1,0) and (2,0) is examined 
in order to to see whether the calculated correlation curve can describe the observational one 
between the upper kHz QPO and HBO found by Psaltis et al. (1999). Two cases of ?7 S = 1.3 
and ?7 S = 1.1 are shown in figure 8 by changing (3 from 1.0 to 20.0. The mass of the central 
source is M = 1.4M and the spin parameter a* = 0. The abscissa is the frequency of the (1,0) 
oscillation, and the ordinate is that of the (2,0) oscillation. The curves are superposed on the 
frequency-frequency diagram plotting observational data points (figure 2.9 in a review paper 
by van der Klis 2004). The blue curve is for ?7 S = 1.3 and the red one is for r] 8 = 1.1 The left-end 
of these curves is for (3 = 20.0 and the right-end is for (3 = 1.0. For comparison, the correlation 
curve between the (1,0) and (2,1) oscillations is also shown by green curve in the case of ?7 S = 1.1 
(in this case the ordinate is the frequency of (2,1) oscillation.) 

We are not interested here in the HBO harmonics and the sub HBO, since their amplitude 
are small and will be subsidiary phenomena. In the case of r] s = 1.3 the correlation curve 
calculated is roughly parallel to the observed kHz QPO - HBO correlation curve, but runs 
much above the observed one on the diagram. As ?7 S decreases the curve moves downward on 
the diagram. The figure shows that in the case of r) s = 1.1 the correlation curve is close to the 
observed sequence of the correlation. We should notice, however, that in the present paper we 
have calculated the frequency correlation of (1,0) and (2,0) oscillations only in a particular case 
of no magnetic fields. If disks have strong magnetic fields, ?7 S required to describe observations 
will not be so small as r] s = 1.1 and also the maximum value of (3 required will become smaller 
than that shown in figure 8 (see next section). 
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Fig. 8. The frequency correlation between (n, n T ) — (1,0) and (2,0) oscillations for two cases of r/ s =1.3 
(second curve) and 1.1 (third curve). The value of (3 is changed from 1.0 to 30.0. The left-end points of 
the curves are for j3 = 30.0, and the right-end points are for j3 = 1.0. The curves are overlapped on the 
plots of observational data which are a part of figure 2.9 by van der Klis (2004). The value of M and a* 
adopted are M = 1.4M and a* = 0. For comparison, the correlation curve between (n,n T ) = (1,0) and 
(2,1) is shown (lowermost curve). The uppermost curve shows, for comparison, the correlation between 
(n, n T ) = (1,0) and (1,1) oscillations (i.e., twin kHz QPOs) in the case of rj s = 1.1. (Color Online) 

6. Discussion 

Trapping of two-armed (m = 2) vertical p-mode oscillations in vertically isothermal disks 
are examined. The disks are assumed to be terminated at a certain height by hot, low-density 
corona. The height of the termination is taken to be a parameter. 

The above issue has been already examined in the previous papers (Kato 2011a,b, 
2012a,b) in the case where the disks have toroidal magnetic fields, introducing the approxima- 
tion that the oscillations are nearly vertical in the sense that the horizontal motions associated 
with the oscillations are small perturbations over the vertical ones. In the case where disks have 
no magnetic fields, however, the trapping can be examined without the above approximation. 
Hence, returning to this simpler case of no magnetic fields, we re-examine in this paper the 
trapping of the vertical p-mode oscillations of m = 2. The purposes of this paper are thus i) to 
examine the case of no magnetic fields more accuratly without the approximation, and ii) to 
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show that even in the case of no magnetic fields the calculated correlation curve can describe 
the observed one, if disk parameters such as f3 = Cg/(Cg)o and t] s = z 8 /H are taken moderately. 

Figure 1 shows that in the case of oscillations of n — 1, the frequencies calculated without 
the above-mentioned approximation are systematically lower compared with those calculated 
with the approximation. Due to this systematical difference, in the case where the approxi- 
mation is not used, the correlation curve between the (1,0) and (1,1) oscillations shifts on the 
frequency - frequency diagram in the left-downward direction almost along the curve derived 
with the approximation (compare figures 5 and 7). We suppose that the same situation occurs 
even in magnetized disks. That is, if we could calculate the correlation curves without the 
approximation in the case of magnetized disks, they will shift in the left-downward direction 
along the curve derived in paper I. This means that the same correlation range on the diagram 
can be derived by smaller values of /3 and c\/c si if frequencies coud be calculated without 
the approximation. In other words, such strong magnetic fields required in paper I may be 
unnecessary, if we could derive and solve wave equations in magnetized disks without the ap- 
proximation. This is in favor of our model, since the magnetic fields required in paper I seem 
to be somewhat stronger than those usually supposed in accretion disks. 

Next, we examine whether the (2,0) oscillation can be regarded as the horizontal branch 
oscillation (HBO). Figure 8 shows that in non- magnetized disks the frequency correlation be- 
tween the (1,0) and (2,0) oscillations is close to the observed correlation between the upper 
kHz QPO and HBO, if ?7 S is taken to be as small as rj B — 1.1. This suggests that the disks of 
horizontal branch phase of Z-sources are strongly diminished in their vertical thickness. 

Concerning the disks of the horizontal branch, we considered in paper I two possibilities: 
The sequence of horizontal branch is i) a sequence of change of toroidal magnetic fields (figure 
9 of paper I) or ii) a sequence of change of disk temperature with moderate or strong toroidal 
magnetic fields (figure 10 of paper I). The first possibility seems to be less promising, since 
the correlation curve calculated in paper I has a sharper gradient on the frequency-frequency 
diagram, compared with the curve of observations (see figure 9 of paper I). This sharp gradi- 
ent comes from the fact that as the toroidal magnetic fields become strong, the frequency of 
the (2,0) oscillation decreases strongly when n = 2§| Hence, the correlation curve calculated 
by changing the strength of toroidal magnetic fields has a sharp gradient on the frequency- 
frequency diagram. This situation will not be much changed, even if we could calculate the 
frequency of trapped (2,0) oscillations in magnetized disks without using the approximation of 
nearly vertical oscillation. 

3 This comes from the following fact that the frequency of the fast mode (the fast mode of the three MHD 
waves) oscillating in the purely vertical direction is given by (u — mil) 2 — [K, hS (c 2 + c\)/ (c 2 + c\ /2) + 
[see equation (47) of Kato (2012a) or equation (27) of Kato (2011a)]. This means that the frequency u> given 
by oj = 2S1 - [K njB (c* + c\)/(c 2 + c\/2) + l] 1/2 £!^ decreases, as (c 2 + c\)/(c 2 + c\/2) increases from 1 to 2 
by increase of c\/c 2 when n = 2 (K n>s > 1 when n = 2). 
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Compared with this, the correlation curve obtained by changing the disk temperature 
can better describe the observational correlation curve. In the limit of no magnetic fields 
the calculated curve is roughly parallel to the observed correlation curve on the frequency 
- frequency diagram and moves downward on the diagram as rj s decreases (see figure 8). If 
toroidal magnetic fields are present (c^/c 2 being taken to be constant during temperature 
change), the correlation curve calculated will move downward furthermore compared with the 
case of no magnetic field, since an increase of c\/c^ will have a similar effect to a decrease of ?7 S 
[see equation uj = 2Q — [K n>s (c* + c\) / (c 2 + c\/2) + lp/^jE This suggests that the observed 
correlation curve will be described by rj B larger than r] s — 1.1, when there are toroidal magnetic 
fields of moderate strength with c^/c 2 = const. The correlation curve shown in figure 10 of 
paper I will be also consistent with this argument, since the correlation curve in the figure will 
shift in the right- downward direction if we could calculate the correlation curve without using 
the approximation of nearly vertical oscillations. 

The above argument suggests that the horizontal branch of Z-sources is a sequence of 
temperature change in disks whose vertical thickness is strongly diminished. The direction of 
temperature increase is from right to left on the horizontal branch, since the direction of decrease 
of observed QPO frequencies is from right to left on the branch. The temperature change will 
be a result of a change of mass accretion rate. Thus, this disk model of the Z-branch is similar 
to and support the disk model proposed by Church et al. (2006) and Bolucinska-Church et al. 
(2010) from spectral analyses, although the origins of QPOs are different. 

One of possible causes why the disks of horizontal branch are vertically thin will be strong 
radiation from the central source resulting from high mass accretion and disk evaporation 
toward the corona, as Church et al. (2006) and Bolucinska-Church et al. (2010) argue. If 
moderately strong toroidal magnetic fields are present in disks, it will be in favor of our model 
of QPOs, although only the limiting case of no magnetic fields is examined in this paper. It 
is noted here that disks with strong magnetic fields (the plasma /3 being less than unity and 
as small as 0.1) are expected at a stage of transition from optically thin flows (ADAFs) to 
optically thick disks by increase of mass accretion rate (Machida et al. 2006, Oda et al. 2007, 
2009, 2010), although the horizontal branch of Z-sources is different from such stage. 

It should be emphasized that in our model of kHz QPOs the correlation curve between 
the (1,0) and (1,1) oscillations moves almost along a common curve on the frequency - frequency 
diagram for changes of disk parameters, as far as the mass and spin of the central source are 
fixed. On the other hand, the correlation curve between the (1,0) and (2,0) oscillations are 
rather sensitive to the disk parameters such as r] s . The former is in favor of our model, since the 
common curve mentioned above can well describe the observed correlation curve for reasonable 
mass and spin parameter of the central source. However, the latter has two sides for and against 

4 A decrease of rj s increases K n , s , while an increase of c\/c^ increases (c|? + c\)/(c^ + c\/2). Hence, both 
have similar effects on frequency w. 
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our model. If our model is correct, it is a good tool to determine the disk structure of horizontal 
branch of Z-sources. It is, however, unclear why disk thickness (and magnetic fields) in real 
disks is really adjusted to values relevant to describe observations. 

In our model the radius of the inner edge of disks is assumed to be unchanged during 
the disk evolution along the sequence. It is unclear at the present stage of our model whether 
time change of the inner edge of disks is required to describe the observed correlation curve by 
our model. 

The observed amplitude of HBOs is much larger than that of kHz QPOs. Whether this 
can be described by our model is a subject to be discussed further. One of possibilities will be 
the difference of width of the trapped regions. In our model, the trapped region of the (1,0) 
oscillation is only in the innermost region, but that of the (2,0) oscillation is wide when rj s is 
small (see figure 3). This difference will be one of causes of amplitude difference between kHz 
QPOs and HBOs. It is noted that the trapped region of the (1,1) oscillation is also wide (see 
figure 2), but in this case the oscillation has one wavelength in the trapped region. Hence, by 
a cancellation by phase mixing the observed amplitude will not become large. It is also noted 
that a wide trapped region implies that our basic approximation that H and (3 are constant in 
the trapped region will not always be relevant. Radial variations of H and /3 should be taken 
into account to do more detailed comparison with observations. 

Appendix A Relation between Trapping Conditions with and without 
Approximation of Nearly Vertical Oscillations 

In a previous paper (Kato 2012a), we have examined the effects of toroidal magnetic fields 
on frequencies of trapped oscillations in truncated disks. In the study we have introduced an 
approximation that the vertical p-mode oscillations consist of nearly vertical motions and thus 
the horizontal motions associated with them can be taken to be small perturbations over the 
vertical ones. Under this approximation, in addition to the WKB approximation, we obtained 
an equation describing the frequencies of trapped oscillations [see equations (39), (41), and (44) 
in Kato (2012a)] 

[(uj-mQ) 2 -^] 1 / 2 e 1 ^ / 3 \ 

r dr = vr n r +- , (15) 



V±H All* V r 4, 

where A rhS is a dimensionless quantity given by equation (37) in Kato (2012a). Furthermore, e 
is a dimensionless quantity defined by [see equations (12) and (17) in Kato (2012a)] 

e(r) = 1 J , L H 2 - K n>s , (16) 

c s -\- c A 

and is assumed to be e < 1 (Kato 2012a). 

Since in the case of no magnetic fields (c\ = 0), A njS is reduced to 

^■n,s = Kn,s + 1 = K n>s , (17) 



equation (fT5l) can be written as 
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rc Uco - mQ) 2 - h 2 } 1 ' 2 r (cu - mQ) 2 - K n S Q 2 1 l ' 2 , 

dr = n [n T + - . (18) 



-"■n,s i4 _|_ 
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This is the equation corresponding to equation (j!4p in this paper, which is 

[(„ - m fl)2 - K 2]l/2 r { _ mn) 2 _ k Q2 ] 1/2 , 3 x 

? dr = vr n r + - . 19 

c s L (a; — ml 2) 2 J V 4/ 

This means that the perturbation method adopted by Kato (2012a) is valid only when in the 
propagation region of oscillations we have 

(cu - mVt) 2 ~ k n>s tt\ (20) 

or 

uo~mVL- Kl /2 Vt ± . (21) 

This is nothing but the condition of e < 1, as is expected. 

In the case of two armed oscillations, m — 2, with n — 1, we have K n s = 1. Hence, this 
condition is realized when w is comparable with Q. In oscillations with n — 2, K n s = 2 in disks 
with ?7 S = oo and K n>s > 2 in the disks terminated in the vertical thickness (see table 1). Hence, 
this approximation is valid when uo is moderately smaller than Q as is shown in figure 1. 
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